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Abstract 


New  concepts  for  communications-based  tactical  decision  aids,  which  can  be  used  for  low-cost 
tactical  UXV  operations,  are  investigated  to  provide  input  for  future  mission  planning  tools  that 
could  be  developed  for  the  Royal  Canadian  Navy.  The  proposed  tactical  decision  aids  provide 
information  on  where  and  when  communications  between  a  UXV  and  its  base  station  will  be 
effective  and  not  effective  as  a  function  of  the  current  and  future  environmental  conditions.  The 
required  signal  propagation  calculations  use  data  from  an  Environment  Canada  numerical  weather 
prediction  model  and  a  terrain  elevation  database.  A  description  is  provided  of  the  approach  used 
to  calculate  the  base  quantities  (signal-to-noise  ratio,  Eb/N0,  bit  error  rate)  that  are  required  to 
form  the  tactical  decisions  aids.  Example  tactical  decision  aids  are  calculated  for  a  specific  case  to 
demonstrate  how  this  information  could  be  used  for  planning  UXV  routes.  The  example  case 
clearly  shows  regions  where  communications  are  effective  and  not  effective  within  a  realistic 
environment.  Experimental  trials  are  required  to  test  the  hypothesis  that  the  proposed  method  is 
capable  of  accurately  now-casting  and  forecasting  the  times  and  places  for  which  a  UXV  will 
have  effective  communication  with  its  base  station. 


Significance  to  defence  and  security 


The  Royal  Canadian  Navy  (RCN)  makes  extensive  use  of  unmanned  aerial  vehicles  (UAV)  on 
ships  at  sea.  Defence  Research  and  Development  Canada  is  also  experimenting  with  UAVs  and 
unmanned  surface  vehicles  (USV)  to  develop  further  concepts  for  their  use  in  RCN  operations, 
such  as  off-board  jamming.  Communication  with  UAVs  or  USVs  (UXVs)  will  be  lost  if  the  UXV 
travels  too  far  from  the  ship  or  if  the  environmental  conditions  are  not  favourable  for  propagation 
of  the  communication  signals.  In  the  case  of  low-cost  tactical  UXVs,  loss  of  communications  will 
prevent  the  operator  from  controlling  the  function  of  the  UXV,  in  which  case  the  UXV  may  enter 
a  holding  pattern  until  communications  are  restored.  Restoration  of  communications  with  the 
UXV  will  require  the  ship’s  crew  to  undertake  additional  measures  that  can  affect  the  overall 
mission  being  conducted  by  the  ship.  This  research  demonstrates  a  method  for  now-casting  and 
forecasting  where  and  when  communications  with  the  UXV  will  be  effective  and  not  effective  as 
a  function  of  current  and  future  environmental  conditions.  A  mission  planning  tool  based  on  this 
method  would  allow  UXV  paths  to  be  planned  to  avoid  areas  where  communications  loss  is 
predicted  and  to  take  advantage  of  environmental  conditions  that  allow  communications  to  be 
extended  to  areas  beyond  the  expected  operating  range  of  the  UXV.  This  should  benefit  RCN 
operations  and  would  also  be  useful  for  planning  experimental  trials  involving  UXVs. 


DRDC-RDDC-2014-R1 1 


Resume 


De  nouveaux  concepts  d’aide  a  la  prise  de  decisions  tactiques  fondees  sur  les  communications, 
pouvant  servir  aux  operations  des  vehicules  sans  pilote  (VSP)  tactiques  economiques,  sont  etudies 
afm  d’adapter  les  outils  de  planification  de  mission  qui  pourraient  etre  com^us  pour  la  Marine 
royale  canadienne.  Les  aides  a  la  prise  de  decisions  tactiques  proposees  foumissent  des 
renseignements  sur  l’endroit  et  le  moment  ou  les  communications  entre  un  vehicule  sans  pilote  et 
sa  station  de  base  fonctionneront  ou  non  de  maniere  efficace  en  fonction  des  conditions 
environnementales  actuelles  et  futures.  Les  calculs  de  propagation  des  signaux  requis  sont 
effectues  au  moyen  de  donnees  provenant  d’un  modele  numerique  de  predictions 
meteorologiques  d’Environnement  Canada  et  d’une  base  de  donnees  d’ elevation  du  terrain. 
L’approche  utilisee  pour  calculer  les  quantites  de  base  (rapport  signal  sur  bruit,  Eb/N0,  taux 
d’erreur  sur  les  bits)  necessaires  pour  concevoir  les  aides  a  la  prise  de  decisions  tactiques  est 
decrite.  Des  aides  de  ce  type  sont  calculees  en  exemple  pour  un  cas  donne  afm  de  demontrer 
comment  ces  renseignements  pourraient  servir  a  planifier  les  routes  des  vehicules  sans  pilote.  Le 
cas  en  exemple  indique  clairement  les  regions  ou  les  communications  sont  efficaces,  et  celles  ou 
elles  ne  le  sont  pas,  dans  un  environnement  realiste.  Des  essais  experimentaux  doivent  etre 
realises  pour  verifier  l’hypothese  selon  laquelle  la  methode  proposee  est  capable  d’etablir,  de 
maniere  precise,  des  evaluations  en  temps  reel  et  des  previsions  du  moment  et  de  l’endroit  ou  les 
communications  entre  un  vehicule  sans  pilote  et  sa  station  de  base  s’effectueront  de  maniere 
efficace. 


Importance  pour  la  defense  et  la  securite 


La  Marine  royale  canadienne  (MRC)  utilise  intensivement  les  vehicules  aeriens  sans 
pilote  (VASP)  a  bord  de  navires  en  mer.  Recherche  et  developpement  pour  la  defense  Canada 
effectue  aussi  des  essais  avec  les  VASP  et  les  vehicules  de  surface  sans  pilote  (VSSP)  pour 
elaborer  des  concepts  plus  pousses  qui  serviraient  aux  operations  de  la  MRC,  comme  le 
brouillage  deporte.  La  communication  avec  un  VASP  ou  un  VSSP  (VSP)  est  perdue  si  le  VSP 
s’eloigne  trop  du  navire  ou  si  les  conditions  environnementales  ne  favorisent  pas  la  propagation 
des  signaux  de  communication.  Dans  le  cas  des  VSP  tactiques  economiques,  la  perte  des 
communications  empeche  l’operateur  de  controler  les  fonctions  du  vehicule  sans  pilote,  auquel 
cas  le  vehicule  peut  effectuer  une  procedure  d’attente  jusqu’a  ce  que  les  communications  soient 
retablies.  Pour  retablir  les  communications  avec  le  vehicule  sans  pilote,  P  equipage  du  navire  doit 
prendre  des  mesures  supplementaires,  lesquelles  peuvent  avoir  une  incidence  sur  l’ensemble  de  la 
mission  du  navire.  L’etude  demontre  une  methode  qui  permet  d’evaluer  la  situation  en  temps  reel 
et  de  prevoir  le  moment  et  l’endroit  ou  les  communications  avec  un  vehicule  sans  pilote 
s’effectueront  ou  non  de  maniere  efficace,  en  fonction  des  conditions  environnementales  actuelles 
et  futures.  Un  outil  de  planification  de  mission  fonde  sur  cette  methode  permettrait  de  planifier  le 
parcours  que  doivent  suivre  les  VSP  afm  d’eviter  les  endroits  ou  Ton  anticipe  la  perte  de 
communications  et  de  profiter  des  conditions  environnementales  qui  favorisent  P extension  des 
communications  a  des  zones  se  situant  au-dela  de  la  portee  operationnelle  prevue  du  vehicule.  11 
profiterait  aux  operations  de  la  MRC  et  servirait  egalement  a  planifier  les  essais  experimentaux 
effectues  sur  les  VSP. 


n 


DRDC-RDDC-2014-R1 1 


Table  of  contents 


Abstract . i 

Significance  to  defence  and  security . i 

Resume . ii 

Importance  pour  la  defense  et  la  securite . ii 

Table  of  contents . iii 

List  of  figures . iv 

List  of  tables . v 

1  Introduction . 1 

2  Overview  of  the  mission  planning  tool  concept . 2 

3  Calculation  of  the  data  fields  (SNR,  Eb/N0,  BER) . 5 

4  Example  case . 9 

5  Summary  and  conclusions . 16 

References . 17 

List  of  acronyms . 18 


iii 


DRDC-RDDC-2014-R1 1 


List  of  figures 


Figure  1 :  Conceptual  map-based  decision  aid  for  UXV  mission  planning.  The  map  centre 

corresponds  to  91.369069°W  longitude  and  74.4378 13°N  latitude . 3 

Figure  2:  SNR  interpolated  as  a  function  of  UXV  range . 4 

Figure  3:  Information  flow  for  SNR,  Eb/N0,  and  BER  calculations . 5 

Figure  4:  Example  GEM  LAM  data  in  HPAC  format . 1 1 

Figure  5:  SNR  that  would  be  present  in  the  UAV  receiver  at  1  km  above  sea  level . 12 

Figure  6:  Map-based  tactical  decision  aid  showing  communications  effectiveness .  14 

Figure  7:  Tactical  decision  aid  showing  communications  effectiveness  along  the  UAV  flight 

path . 15 


IV 


DRDC -RDDC-20 1 4-R 1 1 


List  of  tables 


Table  1:  Base  station  and  signal  parameters . 9 

Table  2:  UAV  and  signal  parameters . 10 

Table  3:  Environment  parameters . 10 


DRDC-RDDC-2014-R.il 


v 


This  page  intentionally  left  blank. 


DRDC-RDDC-2014-R1 1 


1  Introduction 


The  Royal  Canadian  Navy  (RCN)  makes  extensive  use  of  unmanned  aerial  vehicles  (UAV)  on 
ships  at  sea.  Defence  Research  and  Development  Canada  (DRDC)  is  also  experimenting  with 
UAVs  and  unmanned  surface  vehicles  (USV)  to  develop  further  concepts  for  their  use  in  RCN 
operations,  such  as  off-board  jamming.  Communication  with  UAVs  or  USVs  (UXVs)  will  be  lost 
if  the  UXV  travels  too  far  from  the  ship  where  the  strength  of  the  communications  signal  is  too 
low  to  be  distinguished  from  noise.  The  strength  of  the  UXV  communications  signal  can  also  be 
greatly  affected  by  anomalous  propagation  conditions  such  as  radio-frequency  (RF)  ducting, 
which  can  decrease  or  increase  the  strength  of  the  communications  signal.  In  the  case  of  low-cost 
tactical  UXVs,  loss  of  communications  will  prevent  the  operator  from  controlling  the  function  of 
the  UXV,  in  which  case  the  UXV  may  enter  a  holding  pattern  until  communications  are  restored. 
Restoration  of  communications  with  the  UXV  will  require  the  ship’s  crew  to  undertake  additional 
measures  that  can  affect  the  overall  mission  being  conducted  by  the  ship.  The  aim  of  this  research 
is  to  develop  a  concept  for  a  mission  planning  tool  that  can  be  used  to  now-cast  and  forecast 
where  and  when  communications  with  the  UXV  will  be  effective  and  not  effective  as  a  function 
of  current  and  future  environmental  conditions.  Such  a  tool  would  allow  UXV  paths  to  be 
planned  to  avoid  areas  where  communications  loss  is  predicted  and  to  take  advantage  of 
environmental  conditions  that  allow  communications  to  be  extended  to  areas  beyond  the  expected 
operating  range  of  the  UXV.  This  should  benefit  RCN  operations  and  would  also  be  useful  for 
planning  experimental  trials  involving  UXVs.  This  report  provides  a  conceptual  description  of  the 
proposed  tactical  decision  aids  (TDA)  that  would  be  the  basis  of  the  mission  planning  tool  and  a 
description  of  how  the  base  quantities  used  to  form  the  TDAs  are  calculated.  Example  TDAs  are 
generated  for  environmental  conditions  corresponding  to  a  real  data  set. 
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2  Overview  of  the  mission  planning  tool  concept 


A  mission  planning  tool  is  expected  to  use  various  types  of  graphical  displays  (TDAs)  that 
provide  actionable  information  to  the  user  in  an  understandable  way.  Two  types  of  displays  will 
be  discussed  in  this  report.  The  first  is  a  map-based  display  that  presents  the  information  in  the 
context  of  the  local  geography  and  the  second  is  an  “x-y”  plot  format  that  provides  more 
quantitative  information. 

The  forecast  capability  of  the  mission  planning  tool  will  be  provided  through  the  use  of  data  from 
a  re-locatable  mesoscale  numerical  weather  prediction  (NWP)  model.  Such  a  model  provides 
predictions  of  the  thermodynamic  properties  of  the  atmosphere  (e.g.,  pressure,  temperature, 
relative  humidity),  typically  for  a  12  hour  period,  for  a  region  spanning  a  few  hundred  kilometres, 
and  on  a  grid  having  a  horizontal  resolution  on  the  order  of  a  few  kilometers.  The  Environment 
Canada  Global  Environment  Multiscale  Limited  Area  Model  (GEM  LAM)  has  been  used  to 
provide  this  type  of  data  to  DRDC  in  the  past  for  research  being  conducted  within  the  ABCANZ 
(America  Britain  Canada  Australia  New  Zealand)  Focused  Information  Exchange  Group  (FIEG) 
on  infrared,  visual,  and  RF  propagation.  GEM  LAM  could  be  used  to  calculate  a  grid  of 
predictions  centred  in  space  and  time  on  the  operations  area  of  an  RCN  ship.  These  data  could  be 
provided  to  a  ship  at  sea  using  the  Canadian  Forces  Weather  Office  Site  (CFWOS)  exactly  as 
other  meteorological  data  are  currently  provided. 

When  appropriately  combined  with  a  surface  boundary  layer  model  and  a  terrain  elevation 
database,  the  NWP  data  can  be  used  to  calculate  a  field  of  atmospheric  refractivity  in  space 
and  time.  The  refractivity  data,  a  land  cover  database,  and  information  about  the  UXV 
communications  system  can  be  used  to  calculate  the  propagation  of  the  UXV  communications 
signals.  Link  analysis  can  then  be  used  to  determine  the  signal  power  that  would  be  received  by 
the  UXV  as  a  function  of  space  and  time.  If  effective  communication  can  be  characterized  by,  for 
example,  a  particular  signal-to-noise  ratio  being  present  in  the  UXV  communications  receiver, 
then  these  calculations  will  forecast  where  and  when  effective  communications  with  the  UXV 
will  be  possible. 

A  conceptual  example  of  the  map-based  display  is  shown  in  Figure  1.  For  each  forecast  hour,  a 
data  field  would  be  calculated  as  a  function  of  all  three  spatial  dimensions.  The  data  field  could 
be  the  signal-to-noise  ratio  (SNR)  that  would  exist  in  the  UXV  receiver,  or  it  could  be  the  SNR 
that  exists  in  the  base  station  receiver.  Alternatively,  in  the  case  of  digital  communications,  the 
calculated  data  field  could  be  Eb/N0  (ratio  of  energy  per  bit,  Eb,  to  the  spectral  noise  density,  N0) 
or  the  bit  error  rate  (BER).  Figure  1  could  represent  a  horizontal  slice  through  this  field  at  a 
particular  altitude  or  perhaps  an  average  representation  of  the  data  field  for  a  layer  of  altitude.  A 
vertical  range-height  slice  through  the  data  field  along  a  particular  bearing  could  also  be  a  useful 
product.  The  colours  are  used  to  distinguish  where  communications  will  be  effective  and  where 
they  will  not  be.  Green  colour  corresponds  to  areas  where  the  calculated  quantity  (SNR,  Eb/N0, 
or  BER)  exceeds  that  required  for  effective  communications  and  red  colour  corresponds  to  areas 
where  the  calculated  quantity  does  not  exceed  that  required  for  effective  communications.  In  the 
example  shown  in  Figure  1,  the  non-circular  shape  of  the  green  area  would  be  caused  by  effects 
such  as  terrain  shadowing  and  RF  ducting.  A  proposed  UXV  path,  defined  by  a  set  of  waypoints 
is  depicted  with  the  white  line.  The  display  postulates  conditions  such  that  the  range  of  effective 
communications  is  extended  beyond  normal  expectations  in  areas  along  the  coast  of  Devon 
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Island,  Nunavut.  This  could  be  caused,  for  example,  by  dry  air  from  above  the  land  being 
advected  over  the  ocean,  which  could  create  an  elevated  RF  duct.  The  UXV  path  takes  advantage 
of  these  conditions.  However,  the  display  also  shows  that  communications  will  be  lost  at 
waypoint  8  and  that,  therefore,  the  chosen  path  should  be  adjusted  somewhat. 
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Figure  1:  Conceptual  map-based  decision  aid  for  UXV  mission  planning.  The  map  centre 
corresponds  to  91. 369069°W longitude  and  74. 4378 13°N  latitude. 

The  calculated  data  field  could  also  be  interpolated  along  the  UXV  path  and  displayed  in  a  more 
quantitative  fashion  as  a  function  of  variables  such  as  range,  bearing  (azimuth),  altitude, 
longitude,  latitude,  and  time.  For  example,  the  SNR  corresponding  to  the  UXV  path  displayed  in 
Figure  1  is  plotted  as  a  function  of  range  in  Figure  2.  Here  a  value  of  0  dB  has  been  chosen  as  the 
SNR  required  for  effective  communication.  Note  that  the  SNR  values  corresponding  to  the 
waypoints  may  be  calculated  using  NWP  data  from  different  forecast  hours,  depending  upon  the 
time  at  which  the  UXV  reaches  a  given  waypoint. 
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Figure  2:  SNR  interpolated  as  a  function  ofUXV  range. 


If  the  transmission  of  the  GEM  LAM  data  to  a  ship  at  sea  presents  a  problem,  for  example,  due  to 
low  available  data  transfer  rates,  then  the  mission  planning  products  could  be  generated  using  a 
reach-back  process.  This  would  involve  the  ship  transmitting  the  necessary  data  to  shore,  such  as 
ship  location,  proposed  UAV  flight  path,  etc.  The  mission  planning  products  (TDAs)  could  then 
be  calculated  on  shore  more  quickly  using  a  powerful  computing  system,  and  only  the  final 
products  would  be  transmitted  back  to  the  ship.  It  is  expected  that  the  data  size  of  the  final 
mission  planning  products  will  be  much  less  than  the  size  of  the  GEM  LAM  data  grid.  This 
approach  removes  the  requirements  for  computing  power  at  sea  and  for  ship’s  personnel  to 
perform  the  calculations. 


4 


DRDC-RDDC-2014-R1 1 


3  Calculation  of  the  data  fields  (SNR,  Eb/N0,  BER) 


The  information  flow  corresponding  to  the  data  field  calculations  is  depicted  in  Figure  3.  The 
black  boxes  on  the  left  show  required  input  information.  The  specific  information  that  is 
represented  by  each  box  is  described  in  the  parameter  tables  given  for  the  examples  of  Section  4. 
NWP  data  files  may  contain  data  for  more  than  one  forecast  hour.  The  forecast  hour 
corresponding  to  the  calculations  is  determined  based  on  the  proposed  UXV  launch  time,  the 
UXV  speed,  and  the  UXV  location  relative  to  the  base  station. 


Base  Station  Properties 


UXV  Properties 


NWP  Model 
(GEM-LAM) 


Surface  Boundarv  Laver  Model 
(LWKD) 


Terrain  Elevation  Database 
(DTED) 


Land  Cover  Database 


Signal  Properties 


Propagation 
*|  Calculations 
(APM) 


BER 


>  SNR 


Figure  3:  Information  flow  for  SNR,  Eb/Na,  and  BER  calculations. 

The  quantities  Eb/N0  and  BER  are  both  dependent  upon  the  pre-detection  SNR  (see  [1],  p  105). 
The  SNR  in  the  UXV  receiver  is  calculated  as 


SNRU  = 


N„ 


(1) 


where  Pru  is  the  received  power  given  by 


P  = 

rru 


ptb  Gb  flb  Gu  f£u  A2  Ft 


(47t  r)2 


(2) 


and  Ptb  is  the  power  transmitted  by  the  base  station,  Gb  is  the  boresight  gain  of  the  base  station 
antenna,  f£b  is  the  azimuth  pattern  factor  for  the  base  station  antenna,  Gu  is  the  boresight  gain  of 
the  UXV  antenna,  f£u  is  the  azimuth  pattern  factor  for  the  UXV  antenna,  A  is  the  wavelength  of 
the  transmitted  signal,  Fq  is  the  one-way  pattern  propagation  factor,  and  r  is  the  slant  range 
between  the  base  station  antenna  and  the  UXV  antenna. 


DRDC-RDDC-2014-R1 1 


5 


The  noise  power,  Nu,  is  given  by 


Nu  =  kb  T0  Bru  Fnu  (3) 

where  kb  (=  1.38  x  10  23  J/K)  is  the  Boltzmann  constant,  T0  (=  290  K)  is  a  reference  temperature, 
and  Bru  and  Fnu  are  the  bandwidth  and  noise  figure,  respectively,  of  the  UXV  receiver. 

The  pattern  propagation  factor,  Fg,  is  defined  as  the  ratio  of  the  electric  field  that  actually 
impinges  on  a  scatterer  to  the  electric  field  that  would  exist  at  the  scatterer  under  free-space 
conditions  with  the  antenna  positioned  such  that  the  scatterer  is  in  the  maximum-gain  direction  of 
the  antenna  pattern  (see  [2],  p  8-12).  This  factor  can  account  for  the  effects  of  diffraction, 
refraction,  reflection  (multipath  interference),  absorption  by  atmospheric  gases,  surface 
roughness,  and  the  gain  pattern  of  the  antenna.  Use  of  the  pattern  propagation  factor  in  Equation  2 
requires  Fd  to  be  squared  to  relate  the  quantity  to  power. 

In  this  study,  Fg  is  calculated  using  version  2.3.03  of  the  Advanced  Propagation  Model  (APM) 
[3],  The  subscript  8  is  used  to  denote  that  these  calculations  factor  in  the  normalized  polar  angle 
antenna  gain  pattern.  The  APM  calculations  also  include  attenuation  due  to  atmospheric  gases. 
APM  employs  parabolic  equation  methods  [4]  to  calculate  Fg  on  a  range-height  grid  (constant 
azimuth  plane)  as  a  function  of  the  specified  environment.  To  generate  the  three-dimensional  data 
set  of  Fg  that  is  required  for  the  map-based  display,  APM  calculations  are  executed  along  many 
different  azimuths  from  the  ship  to  the  maximum  range  of  interest.  To  generate  the  data  set  of  Fg 
that  is  required  for  the  “x-y”  plot  displays,  APM  calculations  are  executed  along  the  azimuth 
planes  to  a  set  of  points  that  are  evenly  spaced  in  time  along  the  UXV  flight  path. 

The  environment  within  the  range-height  grid  corresponding  to  a  given  path  is  defined  using  a 
number  of  parameters  that  are  inputted  into  the  APM  software.  These  parameters  define  surface 
absolute  humidity,  surface  air  temperature,  surface  wind  speed  as  a  function  of  range,  refractivity 
as  a  function  of  range  and  height,  surface  relative  permittivity  as  a  function  of  range,  surface  ionic 
conductivity  as  a  function  of  range,  and  terrain  height  ( ht )  as  a  function  of  range.  Terrain  height 
is  determined  for  the  examples  of  this  document  by  interpolating  Digital  Terrain  Elevation  Data 
(DTED,  provided  on  compact  disc  by  the  U.S.  National  Imagery  and  Mapping  Agency). 

The  atmospheric  parameters  required  by  APM  are  derived  from  the  GEM  LAM  data.  The 
modified  refractivity  required  by  APM  is  calculated  as  a  function  of  range  and  height  using  the 
Debye  equation  (see  [4],  p.  42): 


M(r,  h ) 


77.6 

Ta (r,  h)  +  273.15. 


P  (r,  h ) 


[4807  es(r,  h)  fr(r,  h)]') 
+  [Ta(r,  h)  +  273.15]  j 


(4) 


where  h  is  height,  Ta,  is  air  temperature  in  degrees  Celsius,  P  is  pressure  in  millibars,  and  fr  is  the 
relative  humidity  expressed  as  a  fraction.  The  saturation  vapour  pressure  over  water  (in  millibars) 
is  given  by  (see  [5],  p  16) 


es  =  6.112  exp 


l  17.67  Ta  \ 
\Ta  +  243.5/ 


(5) 
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The  pressure,  temperature,  and  relative  humidity  data  used  in  Equations  4  and  5  are  calculated  as 
height  profiles  that  are  created  by  merging  (gradient-matching  method)  height  profiles  produced 
by  LWKD  [6]  with  upper  air  data  interpolated  from  the  GEM  LAM  data.  LWKD  is  a  bulk 
aerodynamic  model  that  describes  the  surface  layer  of  the  boundary  layer  using  standard 
meteorological  measurements.  The  boundary  layer  is  the  portion  of  the  atmosphere  that  is 
influenced  by  the  surface  of  the  earth.  The  surface  layer  is  the  lowest  portion  of  the  boundary 
layer  and  it  is  adjacent  to  the  earth’s  surface.  LWKD  requires  input  values  of  water  surface 
temperature,  air  temperature,  air  humidity,  wind  speed,  and  pressure.  All  of  these  inputs  are  also 
interpolated  from  the  GEM  LAM  data. 

The  electrical  characteristics  (relative  permittivity  and  ionic  conductivity)  of  the  surface  over 
which  the  communications  signal  propagates  are  a  function  of  the  specific  terrain  type  (land 
cover).  APM  includes  data  provided  by  the  International  Radio  Consultative  Committee  (CCIR, 
see  Recommendation  527-3  [7])  that  describes  the  relative  permittivity  and  ionic  conductivity  of 
several  general  terrain  types  as  a  function  of  frequency.  For  this  study,  all  land  areas  within  the 
region  of  interest  were  assumed  to  be  “medium  dry  ground”.  A  land  area  is  defined  as  a  region  in 
which  the  terrain  elevation  is  greater  than  zero  with  respect  to  sea  level. 

Depending  upon  the  specifications  of  the  transmitter  and  receiver,  it  is  possible  that  the  limiting 
factor  for  effective  communications  may  be  the  SNR  in  the  base  station  receiver  and  not  the  SNR 
in  the  UXV  receiver.  The  signal-to-noise  ratio  in  the  base  station  receiver,  SNRb,  is  calculated  in 
the  same  way  as  SNRU,  i.e., 


SNRb 


(6) 


where 


Prb  — 


PtuGuf$uGbf£q2F2e 
(47t  r)2 


(7) 


and 


~  k-b  To  Brb  Fnb  (8) 

Here,  Ptu  is  the  power  transmitted  by  the  UXV,  and  Brb  and  Fnb,  are  the  bandwidth  and  noise 
figure,  respectively,  of  the  base  station  receiver.  Assuming  equivalent  antenna  geometry  (e.g., 
omnidirectional  antennas),  the  Helmholtz  reciprocity  theorem  (see  [8],  p  38 1)  dictates  that  the 
value  of  Fq  is  the  same  when  calculated  along  the  path  from  the  base  station  to  the  UXV  as  it  is 
when  calculated  along  the  reverse  path. 
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The  quantity  Eb  / N0  can  be  calculated  from  the  SNR  as 


Eb  B 

-j-  =  SNR  — 
No  fb 


(9) 


where  B  is  the  receiver  bandwidth,  and  fb  is  the  bit  rate  of  the  digital  signal  (see  [1],  p  105).  The 
bit  error  rate,  BER,  can  then  be  determined  from  Eb/N0  depending  upon  the  modulation 
technique  used  to  encode  the  digital  information  onto  the  RF  carrier.  For  example,  for  the  case  of 
quadrature  phase  shift  keying  (QPSK)  modulation  in  additive  white  Gaussian  noise  (AWGN), 
BER  can  be  calculated  as 


1 

BER  =  -erfc 
2 


(10) 


where  erfc  represents  the  complementary  error  function  (see  [9],  p  390).  These  relations  allow  the 
threshold  for  determining  effective  communications  to  be  based  on  the  BER.  For  example,  if 
effective  communications  can  be  achieved  only  if  BER  <  1  X  10-4,  then  assuming  QPSK  and 
AWGN,  and  that  the  bit  rate  equals  the  receiver  bandwidth,  the  required  SNR  is  8.4  dB  or  greater. 
In  this  case  the  map-based  display  (e.g.,  see  Figure  1)  would  be  painted  such  that  green  areas 
correspond  to  SNR  >  8.4  dB  and  red  areas  correspond  to  SNR  <  8.4  dB. 
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4  Example  case 


The  geographic  area  chosen  for  the  example  calculations  is  the  southern  portion  of  the  Delmarva 
Peninsula  and  the  surrounding  waters.  This  choice  was  made  since  GEM  LAM  data  containing 
interesting  propagation  conditions  in  this  area  had  been  previously  provided  to  DRDC  by 
Environment  Canada’s  Recherche  en  Prevision  Numerique  (RPN)  for  studies  being  performed 
within  the  ABCANZ  FIEG. 

It  is  assumed  that  the  UXV  is  a  UAV,  and  that  it  uses  a  radio  system  for  communications  that  is 
identical  to  the  base  station  radio  system  that  is  used  by  its  operator.  The  operator  and  base  station 
are  assumed  to  be  on  a  ship.  The  parameters  associated  with  the  base  station  and  its  transmitted 
signal  are  listed  in  Table  1. 

The  UAV  flight  and  radio  system  parameters  are  listed  in  Table  2.  Waypoint  times  are  shown 
relative  to  the  departure  time  and  are  rounded  to  the  nearest  second. 

The  environment  parameters  associated  with  the  example  calculations  are  listed  in  Table  3.  Based 
on  the  UAV  departure  and  arrival  times,  the  GEM  LAM  data  at  forecast  hour  17:00  UTC  is  the 
most  relevant  for  all  points  along  the  UAV  flight  path.  The  sea  surface  salinity  is  used  in  the 
LWKD  calculations  that  determine  the  thermodynamic  profiles  (temperature,  pressure,  relative 
humidity)  in  the  surface  boundary  layer.  The  temperature,  pressure,  and  relative  humidity  values 
listed  in  Table  3  are  only  used  for  calculating  a  global  value  for  the  gaseous  absorption 
attenuation  rate.  The  temperature,  pressure,  and  relative  humidity  values  that  are  used  to  calculate 
atmospheric  refractivity  are  extracted  from  the  GEM  LAM  data. 


Table  1:  Base  station  and  signal  parameters. 


Description 

Value 

Longitude 

75.3°W 

Latitude 

37.8°N 

Height  above  sea  level 

30  m 

Transmitted  power 

1  W 

Transmitted  frequency 

900  MHz 

Polarization 

Horizontal 

Receiver  bandwidth 

1  MHz 

Receiver  noise  figure 

6dB 

Antenna  type 

Omnidirectional  (f£b  =  1) 

Antenna  gain 

OdB 
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Table  2:  UAV  and  signal  parameters. 


Description 

Value 

Departure  time 

16:31:00  UTC 

Arrival  time 

17:25:46  UTC 

Waypoint  1  (longitude,  latitude,  height,  time) 

75.3°W,  37.8°N,  10m,  Os 

Waypoint  2  (longitude,  latitude,  height,  time) 

75.219543°W,  37.993362°N,  500  m,  227  s 

Waypoint  3  (longitude,  latitude,  height,  time) 

75.219543°W,  38.096401°N,  1000  m,  341  s 

Waypoint  4  (longitude,  latitude,  height,  time) 

75.501770°W,  38.31 1432°N,  1000  m,  685  s 

Waypoint  5  (longitude,  latitude,  height,  time) 

75.492813°W,  38.427906°N,  1000  m,  815  s 

Waypoint  6  (longitude,  latitude,  height,  time) 

75.33 1535°W,  38.414467°N,  1000  m,  956  s 

Waypoint  7  (longitude,  latitude,  height,  time) 

75.341499°W,  37.746975°N,  1000  m,  2099  s 

Waypoint  8  (longitude,  latitude,  height,  time) 

75.278778°W,  37.612579°N,  1000  m,  2258  s 

Waypoint  9  (longitude,  latitude,  height,  time) 

75.4131 17°W,  37.536423°N,  1000  m,  2404  s 

Waypoint  10  (longitude,  latitude,  height,  time) 

75.062752°W,  37.84553 1°N,  500  m,  3072  s 

Waypoint  1 1  (longitude,  latitude,  height,  time) 

75.3°W,  37.8°N,  10  m,  3286  s 

Transmitted  power 

1  W 

Transmitted  frequency 

900  MHz 

Polarization 

Horizontal 

Receiver  bandwidth 

1  MHz 

Receiver  noise  figure 

6dB 

Antenna  type 

Omnidirectional  (f£u  =  1) 

Antenna  gain 

OdB 

Table  3:  Environment  parameters. 


Description 

Value 

Forecast  hour 

17:00  UTC,  May  1,2000 

Sea  surface  salinity 

32  g/kg 

Terrain  type 

“medium  dry  ground”  [7] 

Surface  air  temperature  (attenuation  by  gaseous  absoiption) 

15  °C 

Surface  atmospheric  pressure  (attenuation  by  gaseous  absorption) 

1013.25  mb 

Surface  relative  humidity  (attenuation  by  gaseous  absoiption) 

90% 

10 
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Figure  4  shows  a  portion  of  the  GEM  LAM  data  in  the  Hazard  Prediction  and  Assessment 
Capability  (HPAC)  format.  The  17:00  UTC  forecast  contains  data  for  43584  locations  with 
horizontal  grid  spacing  of  4  km.  The  header  for  each  location  contains  the  date,  time,  longitude, 
latitude,  and  ground  elevation  in  m,  as  well  as  values  for  the  surface  heat  flux  in  W  m'2,  the 
Monin-Obukhov  length  in  m,  the  friction  velocity  in  m  s'1,  the  temperature  scaling  parameter  in 
°C,  the  specific  humidity  scaling  parameter  in  g  g"1,  the  wind  speed  roughness  height  in  m,  the 
temperature  on  the  surface  in  °C,  and  the  relative  humidity  on  the  surface  as  a  percentage.  The 
remaining  data  for  each  location  consists  of  vertical  profiles  of  height  in  m  (Z  values  in  Figure  4), 
wind  direction  in  degrees,  wind  speed  in  m  s'1,  atmospheric  pressure  in  mb,  temperature  in  °C, 
and  relative  humidity  as  a  percentage.  Although  listed  as  0  m,  the  first  height  level  of  each  profile 
actually  corresponds  to  10  m  above  local  ground  level  for  wind  speed  and  direction  values,  and 
1.5  m  above  local  ground  level  for  pressure,  temperature,  and  relative  humidity  values.  The  data 
for  only  two  locations  are  shown  in  the  figure. 


#  creator:  stephane. gaudreault®. ec. gc. ca 

#  DATE:  2010-02-22  14:31:05  UTC 

#  SOURCE:  GEM  LAM  4 KM  &  hpac.tcl  Vl.O 

#  PROFILE  HEIGHT:  ASL 

#  wind  direction:  true 

#  type:  Norwat  forecast 

#  ANALYSIS:  200005011200 

#  START:  200005011700 

#  end:  200005011700 

#  TIME  REFERENCE:  UTC 

#  mode:  hpac  profile  all  Type  1 


#  locations: 

4  3584 

PROFILE 

ID 

YYVMDD 

HOUR 

LAT 

LON 

ELEV 

HFLUX 

MOL 

u* 

T* 

q* 

zO 

Tsurf 

HOURS 

N 

E 

M 

W/M2 

M 

m/s 

OC 

9/9 

m 

OC 

z 

WDIR 

WSPD 

P 

T 

HUMID 

M 

DEG 

M/S 

MB 

C 

% 

1 

000501 

17 

33.6418 

-77. 5962 

0.00 

5. 59e+01 

-2. 78e+02 

1. 64e-01 

1. 36e-02 

-1. 18e-04 

4. 89e-05 

20.31 

0.00 

235.13 

4.95 

1023.14 

20.60 

80.84 

18.88 

235.41 

5.16 

1020.91 

20.49 

77.73 

55.86 

235. 38 

5.44 

1016. 55 

20.16 

77.62 

105.93 

235. 39 

5.76 

1010.68 

19.73 

77.49 

167. 52 

235.64 

6.27 

1003.48 

19.26 

75.28 

242. 55 

235.20 

6.79 

994.77 

18.72 

71.65 

328.67 

233.76 

7.  53 

984 . 84 

18.  39 

60.09 

427.03 

229.29 

7.72 

973. 59 

17.75 

56.18 

537.95 

224. 34 

7.  55 

961.03 

16.96 

54.49 

660.90 

219.89 

7.10 

947.25 

16.08 

52.  31 

794. 39 

216.93 

6.40 

932.46 

15.08 

51.02 

940. 72 

215.69 

5.64 

916.45 

13.92 

51.  55 

1098.42 

215.44 

4.87 

899.43 

12.62 

52.94 

1267.02 

218. 30 

3.91 

881. 50 

11.24 

53.86 

1446.96 

227.61 

3.00 

862.66 

9.71 

55.27 

1637.79 

244.65 

2.  30 

843.00 

8.07 

56.  39 

1839.08 

277.10 

1.85 

822.64 

6.85 

45.71 

2050.93 

290.46 

2.03 

801.67 

6.67 

21.99 

2275.74 

280.75 

2.72 

779.98 

6.82 

11.60 

2511.05 

280.91 

3.  53 

757.90 

5.95 

16.28 

2 

000501 

17 

33.6771 

-77.6044 

0.00 

4. 92e+01 

-7. 04e+02 

1. 62e-01 

1. 61e-02 

-1. 06e-04 

4. 79e-05 

20.14 

0.00 

233. 32 

4.94 

1023.12 

20.48 

82.00 

18.86 

233.72 

5.17 

1020.89 

20.  39 

78.82 

55.84 

233.92 

5.  54 

1016. 53 

20.08 

78.  33 

105.91 

234.21 

5.97 

1010.65 

19.67 

77.40 

167.48 

234.67 

6.53 

1003.46 

19.23 

74.  50 

242. 51 

234.85 

7.06 

994. 75 

18.74 

69.  50 

328.61 

233.78 

7.76 

984 . 82 

18.49 

56.43 

426.99 

229.62 

7.89 

973. 57 

17.85 

52.85 

537.91 

225.13 

7.65 

961.01 

17.05 

51.66 

660.88 

221.26 

7.15 

947.23 

16.15 

50.29 

794. 39 

218.83 

6.46 

932.44 

15.12 

49.88 

940.  72 

217. 55 

5.75 

916.43 

13.95 

50.98 

1098.44 

216.78 

4.99 

899.41 

12.65 

52.47 

1267.06 

219.61 

4.01 

881.48 

11.27 

53.41 

1447.00 

228. 51 

3.14 

862.64 

9.73 

55.07 

1637.85 

244. 57 

2.45 

842.98 

8.08 

56.19 

1839.14 

275.02 

1.99 

822.62 

6.84 

45.99 

2050.97 

289.08 

2.13 

801.65 

6.64 

22.  38 

2275.76 

280. 54 

2.80 

779.97 

6.79 

11.71 

2511.05 

280.78 

3.60 

757.88 

5.94 

16.23 

RHsurf 

% 

100.00 


100. 00 


Figure  4:  Example  GEM  LAM  data  in  HPAC  format. 


To  reduce  execution  time,  the  calculations  are  performed  in  parallel  using  12  CPU  (central 
processing  unit)  cores.  For  the  example  case,  all  calculations  were  completed  in  less  than 
7  minutes  using  version  8.2  of  Interactive  Data  Language  (IDL)  on  an  HP  Z800  Workstation  (2 
Intel  Xeon  X5660  processors  @  2.80  GHz  and  2.79  GHz)  running  the  Windows  7  Enterprise 
64-bit  operating  system.  The  IDL  code  called  version  2.3.03  of  APM  and  version  8.20  of  LWKD. 
SNR  values  are  calculated  on  a  three-dimensional  grid  in  range,  height,  and  azimuth.  Depending 
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upon  the  time  the  calculations  are  performed  relative  to  the  UAV  flight  time,  the  results  can 
represent  a  now-cast  or  a  forecast  of  the  communications  effectiveness  that  the  UAV  will 
experience  along  its  flight. 

The  calculated  SNR  that  would  be  present  in  the  UAV  receiver  at  an  altitude  of  1 000  m  above  sea 
level  is  displayed  in  Figure  5.  The  slice  through  the  SNR  data  field  at  1000  m  was  chosen  for  this 
display  since  the  UAV  spends  most  of  its  times  at  an  altitude  of  1000  m.  In  Figure  5,  the  solid 
black  line  represents  the  coastline,  the  solid  white  line  represents  the  UAV  flight  path,  and  the 
white  X  symbols  represent  the  defined  UAV  waypoints.  The  range  rings  (dashed  black  circular 
lines)  are  shown  at  20  km  intervals.  The  SNR  data  field  represented  by  the  colours  shows  the 
effects  of  multipath  interference  as  ringed  structures  centred  at  the  ship  location.  The  SNR  value 
is  low  in  the  multipath  nulls  and  high  in  the  multipath  peaks.  Figure  5  also  shows  that  the  SNR 
field  varies  with  azimuth.  This  is  a  result  of  non-uniform  atmospheric  refractivity  and  terrain 
elevation. 
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Figure  5:  SNR  that  would  be  present  in  the  UAV  receiver  at  1  km  above  sea  level. 
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A  similar  picture  to  that  shown  in  Figure  5  can  be  created  from  the  SNR  calculated  to  be  present 
in  the  base  station  receiver.  However,  in  this  case  the  colour  at  any  given  location  represents  the 
SNR  value  that  would  exist  in  the  base  station  receiver  if  the  UAV  were  at  the  given  location.  An 
example  figure  for  the  base  station  receiver  SNR  is  not  shown  since  it  would  be  exactly  the  same 
as  that  for  the  UAV  receiver  SNR,  due  to  the  chosen  parameter  set  and  the  reciprocity  associated 
with  the  propagation  calculations. 

Figure  5  could  be  used  as  a  map-based  tactical  decision  aid  or  it  could  be  simplified  further  based 
on  a  threshold  BER  that  is  required  for  effective  communications.  As  discussed  in  Section  3,  for 
the  case  of  QPSK,  AWGN,  a  bit  rate  equal  to  the  receiver  bandwidth,  and  a  requirement  that 
BER  <  1  X  10-4,  the  SNR  required  for  effective  communications  is  8.4  dB  or  greater.  Figure  6 
shows  the  SNR  data  of  Figure  5  contoured  only  at  8.4  dB.  That  is,  values  greater  than  8.4  dB  are 
shown  as  green  and  values  less  than  8.4  dB  are  shown  as  red.  The  red  and  green  colours  are 
displayed  as  transparent  so  that  the  local  geography  can  also  be  seen.  The  display  shows  that 
effective  communications  are  easier  to  achieve  in  areas  that  are  north-west  and  south-west  of  the 
ship.  The  display  also  suggests  that  communications  may  be  disrupted  at  a  few  points  along  the 
UAV  flight  path  due  to  multipath  interference  and  that  communications  will  be  lost  as  the  UAV 
travels  to  it  farthest  points  from  the  ship  in  the  eastern  direction.  This  type  of  map-based  tactical 
decision  aid  can  be  used  quickly  and  easily  to  determine  that  the  UAV  path  should  be  altered  to 
avoid  loss  of  communications,  and  to  determine  which  areas  the  UAV  flight  path  should  be 
limited  to. 

Within  a  mission  planning  tool,  more  quantitative  information  could  be  added  to  the  tactical 
decision  aid  in  Figure  6  by  creating  a  mouse-over  feature  that  displays  specific  numbers  such  as 
SNR  or  UAV  flight  path  parameters  when  the  user  moves  the  mouse  over  a  particular  location. 
Alternatively,  more  quantitative  information  could  be  provided  through  an  “x-y”  plot  type 
display.  Figure  7  shows  the  SNR  for  both  the  UAV  and  base  station  receivers  as  a  function  of 
time  along  with  the  UAV  flight  path  parameters.  Note  that  the  SNR  for  both  receivers  is  identical 
in  this  case.  Hence,  the  dark  blue  and  light  blue  lines  in  the  SNR  plot  are  coincident.  For  this 
display,  the  SNR  data  does  not  correspond  to  a  specific  altitude  or  azimuth.  Rather  the  data 
correspond  exactly  to  the  UAV  position  for  all  plotted  data  points.  This  type  of  tactical  decision 
aid  will  allow  the  user  to  determine  at  what  times  communications  will  be  effective  and 
ineffective,  and,  for  those  times,  to  read  the  SNR  value  off  the  plot  along  with  the  corresponding 
longitude,  latitude,  range,  azimuth,  and  height  of  the  UAV. 
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Figure  6:  Map-based  tactical  decision  aid  showing  communications  effectiveness. 
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Figure  7:  Tactical  decision  aid  showing  communications  effectiveness  along  the  UAV flight  path. 
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5  Summary  and  conclusions 


This  research  has  proposed  concepts  for  TDAs  that  could  be  used  in  a  mission  planning  tool  to 
now-cast  and  forecast  the  effectiveness  of  communications  between  a  ship  and  a  UXV  as  a 
function  of  environmental  conditions  and  the  planned  UXV  path.  The  methods  used  to  calculate 
the  TDAs  were  described  and  example  TDAs  were  generated  for  real  environmental  conditions. 
The  example  case  showed  the  variation  in  communications  effectiveness  that  can  be  caused  by 
multipath  interference  and  non-uniform  atmospheric  refractivity.  In  addition,  the  example  TDAs 
clearly  indicated  at  which  portion  of  the  UXV  path  communications  would  be  affected  and  how 
the  path  could  be  altered  to  improve  communications  effectiveness. 

The  main  uniqueness  of  this  work  is  the  use  of  re-locatable  mesoscale  NWP  model  data  to 
determine  a  four-dimensional  (space  and  time)  field  of  atmospheric  refractivity  that  can  be  used 
to  calculate  the  propagation  environment  for  the  communications  signals.  The  example  case 
showed  variation  in  communications  effectiveness  as  a  function  of  both  range  and  azimuth.  If  the 
propagation  calculations  were  performed  using  an  atmospheric  refractivity  field  that  was 
horizontally  extrapolated  from  a  single  refractivity  profile,  such  as  that  determined  from  a 
radiosonde,  then  these  range  and  azimuth  variations  in  communication  effectiveness  would  not  be 
observed,  even  though  they  may  actually  exist  in  the  current  or  future  environment.  This  issue 
would  be  most  important  in  a  littoral  environment  where  a  single  refractivity  profile  will  not 
accurately  represent  the  environment  within  the  full  potential  operating  area  of  a  UAV. 

Only  two  example  TDAs  have  been  constructed  in  this  report.  However,  once  the  propagation 
calculations  have  been  performed,  many  different  measures  of  effectiveness  beyond  SNR,  Eb/N0, 
and  BER  could  be  calculated,  and  many  different  types  of  TDA  products  could  be  used  to  relay 
this  information  to  an  operator.  Further  TDAs  along  these  lines,  as  well  as  a  productized  mission 
planning  tool,  should  be  developed  in  consultation  with  the  end  user. 

Initial  plans  for  this  research  included  experimental  trials  to  test  the  hypothesis  that  the  use  of 
NWP  data  in  the  manner  described  can  produce  accurate  now-casts  and  forecasts  of 
communications  effectiveness.  These  trials  would  involve  laboratory  measurements  to  determine 
the  appropriate  threshold  (e.g.,  SNR  level)  that  defines  effective  communications  for  the  UXV 
system  under  investigation  and  field  tests  to  determine  if  communications  are  maintained  and  lost 
as  predicted  (where  and  when)  by  the  generated  TDAs.  However,  the  sponsoring  project  was 
ended  prematurely  due  to  funding  restrictions.  Consequently,  this  research  will  conclude  with  this 
report  unless  further  interest  arises  from  potential  users. 
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DTED 

Digital  Terrain  Elevation  Data 
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Focused  Information  Exchange  Group 

GEM  LAM 

Global  Environment  Multiscale  Limited  Area  Model 

HP  AC 

Hazard  Prediction  and  Assessment  Capability 

IDL 

Interactive  Data  Language 

LWKD 

Luc  (Wavy)  Walmsley  KEL  DREV 

NWP 

Numerical  Weather  Prediction 

QPSK 

Quadrature  Phase  Shift  Keying 

RCN 

Royal  Canadian  Navy 

RPN 

Recherche  en  Prevision  Numerique 

SNR 

Signal  to  Noise  Ratio 

SR 

Scientific  Report 

TDA 

Tactical  Decision  Aid 

UAV 

Unmanned  Aerial  Vehicle 

USV 

Unmanned  Surface  Vehicle 

UTC 

Universal  Time  Constant 

UXV 

Unmanned  Vehicle  of  general  type 
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